Disentangling the dynamical origin of Pn Nucleon Resonances 
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We show that two almost degenerate poles near the ttA threshold and the next higher mass pole 
in the Pn partial wave of nN scattering evolve from a single bare state through its coupling with 
ttN, rjN and ttttA'' reaction channels. This finding provides new information on understanding the 
dynamical origins of the Roper iV*(1440) and A'^*(1710) resonances listed by Particle Data Group. 
Our results for the resonance poles in other ttN partial waves are also presented. 
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The excited nucleon states are unstable and couple 
strongly to the meson-baryon continuum states to form 
resonances in nN and reactions. Therefore, the ex- 
traction of nucleon resonances (called collectively as N*) 
from data has been a well recognized important task in 
advancing our understanding of strong interactions. The 
N* parameters listed and periodically updated by Parti- 
cle Data Group [H (PDG) are commonly used in testing 
hadron structure calculations using QCD-based hadron 
models 0-11] and Lattice QCD [5, 6]. 

It is well known that resonances locate on the un- 
physical sheets of the complex energy plane and thus 
their properties can only be extracted from the empir- 
ical partial-wave amplitudes (PWA) by analytic contin- 
uation. In extracting resonances from nN data up to 
invariant mass W=2 GeV we face a multi-channel com- 
plication, namely that a resonance may appear as a pole 
on more than one of the unphysical Riemann sheets, as 
investigated previously by Eden and Taylor , Kato [1] , 
and Morgan and Pennington It is custom to name 
the pole which is closest to physical region as the res- 
onance pole, and others as shadow poles. In general, 
the observables are mainly determined by the resonance 
poles. However, under certain circumstances a shadow 
pole could lie close to the threshold of one of the chan- 
nels and could therefore affect the physical observables, 
as discussed in Refs. 0, Hi- A theoretical understanding 
of the dynamical origins of these poles and their inter- 
relations is needed to interpret the resonance parameters. 
In this letter, we report a progress in this direction for 
the N* in the Pu partial wave of nN scattering. Our 
results for other partial waves will also be presented. 

The determination of resonance poles in the Pu partial 
wave has been difficult since the discovery [10| of the 
Roper, A^* (1440), resonance in 1964. It was first found by 
Arndt, Ford and Roper |ll| that this partial wave has two 
almost degenerate poles near the ttA threshold. This was 
confirmed and investigated in more detail by Cutkosky 
and Wang [l^l. This two pole structure has also been 



obtained in the recent analysis by the GWU/VPI [13| 
and Jiilich [13] groups. In this letter, we demonstrate 
that these two poles near the ttA threshold 1360 MeV) 
and a pole at about 1800 MeV correspond to a single bare 
state within a dynamical coupled-channels model (JLMS) 
developed in Ref. .15.] . Thus they have the resonance 
pole-shadow pole relation as discussed in Refs. 0-(9] . Our 
result suggest that the N*{U40) and Af*(1710) listed 
by PDG originate from the same excited nucleon state 
modeled as a bare particle within the JLMS model. 

The JLMS model is defined within a Hamiltonian for- 
mulation [16] of multi-channels reactions. It describes 
meson-baryon (MB) reactions involving the following 
channels: nN, r]N, and nnN which has ttA, pN, and 
ctA^ resonant components. The excitation of the internal 
structure of a baryon (B) by a meson (M) to a bare N* 
state is modeled by a vertex interaction Pmb-m-at*- The 
Hamiltonian also has energy independent interactions 
vmb,m'B' which describe the meson-exchange mecha- 
nisms deduced from phenomenological Lagrangians. Nu- 
cleon resonances can be due to the MB N* M'B' 
transitions induced by the vertex interaction Tmb-it^n* 
in this formulation. But they can also be due to the at- 
tractive forces of vmb,m'B' and channel coupling effects. 
For investigating the N* structure, the second type of 
resonances, called molecular-type resonances in the lit- 
erature, must also be identified in the analysis. For the 
same consideration, the parameterization of V]\ib,M'B' , 
in particular their phenomenological form factors, must 
be carefully constrained by the data. This had been 
achieved by performing rather complex y^-fits to the 
nN scattering data, as detailed in Ref. [15|. Briefly, the 
JLMS model is able to describe the data of nN elastic 
scattering up to invariant mass W = 2 GeV. The result- 
ing nN scattering amplitudes and total cross sections are 
in good agreement with those from SAID Further- 



more, the predicted 2n production cross sections [17[ are 
in good agreement with the available data. 

Within the JLMS model, it is convenient to cast the 
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partial-wave amplitude of the M{k) + B{-k) M'{k') + 
B'{—k') reaction into the following form (suppressing the 
angular momentum and isospin indices): 

TMB.M'B'ik,k' ,E) ^ tMBM'B'{k,k' ,E) 

+ i%iB,M' B'{k,k , E) , (1) 

where the first term (called meson-exchange amplitude 
from now on) is defined by 

tMB,M'B'{k,k' ,E) = VMB.M'B'ik,k') (2) 

+ ^ / dq VMB,M"B"{k,q) 

M"B" ■^C'a/"b" 

X Gm"B" (?, E)tM"B",M'B' {q, k'E) , 

where Cmb is the integration contour in the complex— g 
plane used for channel MB. The term associated with 
the bare N* states in Eq. ([T]) is 



tMB,M'B>{k, k', E) 



E 

N'.N* 



fMB^N'{k,E)[D(E)] 



X VN'^M'B'{k' ,E), 



(3) 



where TN*^M'B'ik,E) is the dressed vertex func- 
tion which is calculated [15] from the bare vertex 
TN*-^M'B'{k) and convolutions over the meson-exchange 
amplitudes tMB,M'B'{k, k' , E). The inverse of the prop- 
agator of dressed A^* states in Eq. ^ is 



\D-'(E)\^.J = iE-M%,)S.^,~[MiE)]. 



(4) 



where Mjy. is the bare mass of the z-th N* state, and 
the N* self-energy is defined by 



[MiE)], 



E 



q dqTNj^MB{q,E) 
X Gmb {q, E) Tmb^n* (g, E) . 



(5) 



Defining Ea{k) — [rn^ + fc^]^/^ with TOq being the mass 
of particle a, the meson-baryon propagators in the above 
equations are: GMB{k,E) = 1/[E — EM{k) - Esik) + ie] 
for the stable irN and rjN channels, and GMB{k,E) = 
l/[E - Euik) - Eeik) ~ T.MB{k,E)] for the unstable 
ttA, pN , and aN channels. The self energy T,MB{k,E) 
is calculated from a vertex function defining the decay 
of the considered unstable particle in the presence of a 
spectator tt or N with momentum k. For example, we 
have for the ttA state. 



s.A(fc,£;) 



rriA 



q dq 



M^N{q) 



EA{k)Jc,'' ^[Af2^(g)-fP]V2 

\fA^^N{q)f 

E-E,{k)-[M^^iq) + k^]^/^ + ie 

where MT,iq{q) = ^^^(g) -I- Eiq{q) and fA^wNiq) defines 
the decay of the A ttN in the rest frame of A, C3 is 



the corresponding integration contour in the complex— g 
plane. The self-energies for pN and aN channels are 
similar. 

To search for resonance poles, we need to choose the 
contours Gmb and C3 appropriately to solve Eqs. ([S])-® 
for E on the various possible unphysical sheets of the Rie- 
mann surface. This requires careful examinations of the 
locations of the on-shell momentum of each propagator 
GMB{k,E) and the ttttN cut in the self energies, such 
as T,TrA{k,E) of Eq. ([6]), of the unstable particle chan- 
nels. Furthermore, we need to account for the singulari- 
ties of vmb,m' b' {k, k') of Eq. ([2]) on the chosen contours. 
Our method was tested 18| within several exactly solv- 
able models. Like all previous works [l^, we only 
look for poles which are close to the physical region and 
have effects on -kN scattering observables. All of these 
poles are on the unphysical sheet of the ttN channel, but 
could be on either unphysical (u) or physical {p) sheets 
of other channels considered in this analysis. We will in- 
dicate the sheets where the identified poles are located 
by {st,n, Sr^N, Sr,-^N, SttA, SpN, Sctn), where smb and s^^at 
can be u or p or — denoting no coupling to this channel. 

Eq. ([1]) indicates that if no pole is found in the first 
term t7rAf,7rAr(fc, k', E), then the poles of the total ampli- 
tude can be found from the second term t^^ T^N^k, k', E). 
But if tTrN,TrNik,k' , E) has a pole, we need to check 
whether it will be canceled by the second term, as demon- 
strated in Ref. [31 ■ Thus our procedure is to first use the 
standard method to determine whether t^7v,7r7v(fc, fc', 
has poles by examining the determinant of [1 — vG]~^ 
of Eq. dSJ. It turns out that we don't find any pole 
from these meson-exchange amplitudes. Thus there is no 
molecular-type nucleon resonance within JLMS model. 

We thus can search for poles of the total amplitudes 
from finding the zeros of the determinant of D~^{E) de- 
fined by Eq. ([4]). Here we use the well-established New- 
ton iteration method. We have performed searches in the 
(m^ + mN) < Re{E) < 2000 MeV and -lYn{E) < 250 
MeV region within which PDG's 3- and 4-stars reso- 
nances are listed. Poles with very large widths are more 
difficult to locate precisely with our numerical methods 
and hence will not be discussed here. 

We now focus on our results in Pn partial wave. We 
find two poles near the PDG value (Re Mn, — Im Mfi) = 
(1350-1380,80-110) of the Roper, TV* (1440), resonance. 



TABLE I: Pii resonance pole positions Mr [listed as 
(Re Mu, — Im Mr)] extracted from four different approaches 
are compared. 



Analysis 




Pll poles 


(MeV) 




JLMS fl5] 


(1357, 


76) 


(1364, 


105) 


CMB [12J 


(1370, 


114) 


(1360, 


120) 


GWU/VPI [13| 


(1359, 


82) 


(1388, 


83) 


Jiilich [14] 


(1387, 


74) 


(1387, 


71) 



3 



- TiN.TITC N 





-100 



> 

01 



^ -200 
E 



-300 




6(1364,-105) 



C(1 820,-248) ' 



aN .-- 



1400 1600 1800 

Re (E) (MeV) 



-ImE (MeV 
100 




1300 



1700 



Re E (MeV) 



1900 



TABLE II: The resonance pole positions Mr [listed as 
(Re Mr, — Im Mr)] extracted from the JLMS model in the 
different unphysical sheets are compared with the values of 
3- and 4-stars nucleon resonances listed in the PDG [J. 
The notation indicating their locations on the Riemann sur- 
face are explained in the text. " — " for P33(1600), P13 and 
P31 indicates that no resonance pole has been found in the 
considered complex energy region, Re(£') < 2000 MeV and 
-Im(£;) < 250 MeV. All masses are in MeV. 
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FIG. 1: (above) Trajectories of the evolution of Pn resonance 
poles A (1357,76), B (1364,105), and C (1820,248) from a bare 
A'^* with 1763 MeV, as the couplings of the bare A''* with the 
meson-baryon reaction channels are varied from zero to the 
full strengths of the JLMS model. See text for detailed expla- 
nations. Brunch cuts for all channels are denoted as dashed 
lines. The branch points, -Bb.p., for unstable channels are 
determined by E^.p. — EM(k) — EB(k) — Ems (fe, -Eb.p.) = 
of the their propagators (described in the text) evalu- 
ated at the spectator momentum fc=0. With the param- 
eters [Tl] used in JLMS model, we find that Eb.p. (MeV) 
= (1365.40, -32.46), (1704.08, -74.98), (1907.57, -323.62) for 
ttA, pN, and aN, respectively, (below) 3-Dimensional depic- 
tion of the behavior of \det[D(E)]f of the Pn A* propagator 
(in arbitrary units) as a function of complex-B. 



This finding is consistent with the results from the anal- 
ysis by Cutkosky and Wang [ll] (CMB), GWU/VPI [Ts^ 
and Jiilich [1^] groups, as seen in Tab. HI In our analysis, 
we find that they are on different sheets: (1357,76) and 
(1364,105) are on the un-physical and physical sheet of 
the ttA channel, respectively. 

We also find one higher mass pole at (1820, 248) in 
Pii partial wave, which is close to the iV*(1710) state 
listed by PDG. Within the JLMS model, we find that 
this pole and the two poles listed in table HI] are related 
to one of the two bare states needed to obtain a good 
fit to the Pii amplitude up to = 2 GeV, see p^ . 



To see how these poles evolve dynamically through their 
coupling with reaction channels, we trace the zeros of 
det[D-^(E)] = det[E - M%, - Y.MByMBMMB{E)] in 
the region < umb < 1, where Mmb{E) is the con- 
tribution of channel MB to the self energy defined by 
Eq. ([5]) . Each um b is varied independently to find contin- 
uous evolution paths through the various Riemann sheets 
on which our analytic continuation method is valid. 

We find that the three poles listed in Table Uare asso- 
ciated to the bare state at 1736 MeV as shown in Fig. [T] 
The solid blue curve shows the evolution of this bare 
state to the position at C(1820, 248) on the unphysical 
sheet of the ttA and rjN channels. The poles A(1357, 76) 
and B(1364,105) evolve from the same bare state on the 
physical sheet of the rjN channel. The dashed red curve 
indicates how the bare state evolves through varying all 
coupling strengths except keeping y^rA = 0, to about 
Re(Affl) - 1400 MeV. By further varying y^rA to 1 of the 
full JLMS model, it then splits into two trajectories; one 
moves to pole A(1357,76) on the unphysical sheet and 
the other to B(1364, 105) on the physical sheet of ttA 
channel. Fig. [T] clearly shows how the coupled-channels 
effects induces multi-poles from a single bare state. The 
evolution of the second bare state at 2037 MeV into 
a resonance at W > 2 GeV can be similarly investigated, 
but will not be discussed here. 

To explore this interesting result further and to ex- 
amine the stability of the determined three Pn poles, 
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we have performed several refits of the Pu amphtudes 
within the JLMS model. We are able to get new fits by 
varying solely the parameters associated with the bare 
N* state at 1763 MeV while keeping its bare mass value 
varied within the range 1763 ± 100 MeV. The quality of 
these fits are comparable to that of the original JLMS 
model. The above described features remain unchanged: 
we find in all refitted results two poles close to the ttA 
threshold, within 1 MeV of the positions reported in Ta- 
ble II. The third higher mass pole is also found but its 
position varies up to 30 MeV from the value given in Ta- 
ble II. The trajectories similar to that shown in Fig.[l]are 
also obtained. This is the extent to which the stability of 
the resonance pole-shadow pole relation among the three 
Pii poles we can establish here. A more detailed analysis 
of the model dependence of our results would involve ex- 
tensive refits by varying the parameters associated with 
both the meson-exchange interaction vmb.M' B' and bare 
A''* states in all partial waves and can not be addressed 
here. 

To further compare our Pu poles with the A* (1440) 
and A^*(1710) listed by PDG, we have applied the 
method explained in Ref. [2(i] to extract the residues 
F ~ Re^'^ which is related to S-matrix by S{E) ^■ 
1 -I- 2iF/{E - Mr) a,s E ^ Mr. We obtain 
(i?[MeV],0[degrees]) = (36, -111), (64, -99), (20, -168) 
for the Pu poles at (Re M^, -Im Mr) = (1357,76), 
(1364, 105), and (1820, 248), respectively. The branching 
ratio of the A^* decay into ttN channel can then be esti- 
mated by evaluating rje ^ ^/(— Iiii(Mfl). Our results for 
the Pii poles at (1357, 76) and (1364, 105) are 49% and 
61%, respectively. These values are close to 60 — 70% 
of the N*{1U0) listed by PDG. Our result for the pole 
at (1820, 248) is 8% which is also close to 10 - 20% of 
A'* (1710). We thus have firmer evidence showing that 
these two A^* states listed by PDG do evolve from the 
same bare state through its coupling with ttN, r/N, and 
ttttN reaction channels. 

Let us now turn to other partial waves. In Table |lll 
the extracted resonance poles positions {Mr) are com- 
pared with the bare A^* masses (M^.) of the JLMS 
model and the 3- and 4-star values listed by PDG [H. 
With the exception of the P33(1600), P13 and P31 cases, 
all pole positions listed by the PDG are consistent with 
our results. One possible reason for not finding these 
poles is that their imaginary part may be beyond the 
— Im(Mi{) < 250 MeV region where our analytic contin- 
uation method is accurate and is covered in our searches. 
Another possibility is that these resonances, if indeed ex- 
ist, are perhaps due to the mechanisms which are beyond 
the JLMS model, but are particularly sensitive to these 
partial waves. On the other hand, the possibility that 
these resonances do not exist can not be excluded since 
the ttA' data are not complete and all partial wave analy- 
ses involve unavoidable theoretical assumptions. For the 
P35 partial wave, we have also analyzed the evolution tra- 



jectories and found that the two poles listed in Table [TT] 
correspond to the same bare state at 2162 MeV. 

In summary, we have applied an analytic continuation 
method [isj to extract nucleon resonances from a dy- 
namical coupled-channels model within which the bare 
A^* states were determined from fitting the ttN scatter- 
ing data wp to W — 2 GeV [13] . Compared with all 
previous analysis, the new aspect of this work is to study 
the evolution of resonance pole parameters as a function 
of the coupling to continuum meson-baryon channels. 
Our most important finding is that the two lowest Pu 
nucleon resonances, the Roper A^*(1440) and A^*(1710), 
originate from a single bare state. Our finding has an 
important implication in understanding how nucleon res- 
onances arise in QCD. It implies that in some limits in 
which the coupling to the continuum is not fully imple- 
mented, for example large Nc QCD or quenched lattice 
QCD, there could be fewer nucleon resonances. Another 
possible implication is that the bare A^* states, not the 
resonance poles, determined within our model could cor- 
respond to hadron structure calculations which exclude 
the coupling with meson-baryon continuum. Further in- 
vestigations of these possibilities as well as related theo- 
retical questions are needed to open a new direction to- 
wards understanding nucleon resonances and their con- 
nection to QCD. Finally, we mention that our results 
have confirmed most of the 3- and 4-stars nucleon reso- 
nance poles listed by PDG but found no evidence of two 
four star resonances, Pi3(1720), P3i(1910), and one three 
star one, P33(1600). 
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